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ABSTRACT 


Theoretical 0 IV electron density sensitive emission line ratios, determined 
using electron impact excitation rates calculated with the R-matrix code, are 
presented for Ri = 1(1407.4 A)/I(1401.2 A), R2 = 1(1404.8 A)/I(140 1.2 A), R 3 
= 1(1399.8 A)/I(1401.2 A) and R 4 = 1(1397.2 A)/I(1401.2 A). The observed 
values of Ri - R 4 , measured from high resolution spectra obtained with the 
International Ultraviolet Explorer (IUE) satellite, lead to electron densities 
that are compatible, and which are also in good agreement with those deduced 
from line ratios in other species. This provides observational support for the 
accuracy of the atomic data adopted in the present calculations. 


Subject headings: atomic data — binaries: symbiotic — planetary nebulae: 
general — ultraviolet: stars 
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1. INTRODUCTION 

The O IV 2s*2p 2 P - 2s2p 2 4 P intercombination lines near 1400 A have 
been observed in the spectra of astrophysical sources, including planetary neb- 
ulae (Hayes k Nussbaumer 1983), symbiotic stars (Feibelman 1982), and the 
solar transition region (Sandlin et al. 1986). Flower k Nussbaumer (1975) 
first pointed out the diagnostic applications of emission line ratios involving 
these transitions, and presented relative line strengths determined using elec- 
tron impact excitation rates calculated in the distorted-wave approximation 
(Eissner k Seaton 1972). Hayes k Nussbaumer (1983) subsequently rederived 
these ratios using more accurate electron rates deduced with the IMP ACT 
computer code of Crees, Seaton, & Wilson (1978). 

Recently Blum k Pradhan ( 1992) have calculated excitation rates for tran- 
sitions in 0 IV using the R-matrix method as adapted for the Opacity Project 
(Seaton 1987; Berrington et al. 1987), and found results somewhat different 
from those of Hayes k Nussbaumer (1983). For example, at T e = 5000 K, 
the Blum k Pradhan rate for the 2s 2 2p 2 P 1/2 - 2s2p 2 4 P 1/2 transition is ~15% 
larger than that of Hayes k Nussbaumer, while for the 2 sip 1 4 Pi/ 2 - 2 sip 1 
4 P 5/2 transition the R-matrix data are ~20% larger. In this paper we use the 
Blum k Pradhan results to derive emission line ratios applicable to gaseous 
nebulae, and compare these with observations obtained with the International 
Ultraviolet Explorer (IUE) satellite. 



2. THEORETICAL RATIOS 


The model ion for 0 IV consisted of the 15 energetically lowest fine- 
structure levels, namely 2^2? 2 Pi/ 2 , 3 / 2 ; 2s2 p 2 4 Pi/ 2 ,3/2,s/2, 2 E 3 / 2t5 / 2 , 2 Si/ 2 , 
2 Pi/ 2 , 3 / 2 l 2p 3 4 S 3/2 , 2 D 3/2 , 5/2 and 2 P 1/2 , 3 / 2 - Energies for all these levels were 

obtained from Fawcett (1975). 

Electron impact excitation rates for transitions in 0 IV were taken from 
Blum k Pradhan (1992), while for Einstein A-coefEcients the calculations of 
Nussbaumer k Storey (1982) and Dankworth k Trefftz (1977) were adopted. 
As noted by, for example Seaton (1964), excitation by protons may be im- 
portant for fine-structure transitions. However Flower k Nussbaumer (1975) 
found proton rates for transitions within 2^2? 2 P and 2s2p 2 4 P to be typically 
a factor of 100 smaller than the corresponding electron rates at the tempera- 
tures considered here. Hence these rates should have a negligible effect on the 
theoretical line intensities for 0 IV. 

Using the atomic data discussed above in conjunction with the statistical 
equilibrium code of Dufton (1977), relative 0 IV level populations and hence 
emission line strengths were calculated for a range of electron temperatures 
and densities. The following assumptions were made in the calculations, (i) 
that photoexcitation and de-excitation rates are negligible in comparison with 
the corresponding collisional rates, (ii) that ionization to and recombination 
from other ionic levels is slow compared with bound-bound rates, and (iii) that 
all transitions are optically thin. Further details of the procedures involved 
may be found in Dufton (1977) and Dufton et al. (1978). 

In Figures 1 and 2 we plot the theoretical ratios 

Rj = I(2s 2 2p 2 P 3/2 - 2s2p 2 4 P 1/2 )/I(2 s 2 2 p 2 P 3/2 - 2s2p 2 4 P s/2 ) 
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= 1(1407.4 A)/I(1401.2 A), 

R 2 = 1(2^2? 2 P 3/2 - 2s2p 2 4 P 3/2 )/I(25 2 2i7 2 P 3/2 - 2s2p 2 4 ? s/2 ) 

= 1(1404.8 A)/I(1401.2 A) 

as a function of electron density for a range of electron temperatures, T e = 
5000 - 40000 K. An inspection of the figures reveals that the temperature 
dependence of the ratios is small, but that they vary with density for N e > 
10 2 cm -3 . For example, at T e = 10000 K, R x and R 2 change by factors of 2.8 
and 2.1, respectively, between N e = 10 2 and 10 s cm -3 . However the variation 
in the ratios between T e = 10000 and 20000 K is only 11% (Ri) and 9% (R 2 ) 
at N e = 10 2 cm -3 , which decreases to 2% (Rj) and 1% (R 2 ) at N e = 10 6 cm -3 . 
We note that the ratios 

R 3 = I(2s 2 2p 2 P 1/2 - 2s2p 2 4 P 1/2 )/I(2 s 2 2 p 2 P 3/2 - 2s2p 2 4 P 5/2 ) 

= 1(1399.8 A)/I(1401.2 A), 

R 4 = 1(2^2? 2 P 1/2 - 2s2p 2 4 P 3/2 )/I(2s 2 2p 2 P 3/2 - 2s2p 2 4 P 5/2 ) 

= 1(1397.2 A)/I(1401.2 A) 

have the same density dependence due to common upper levels as Ri and R 2 , 
respectively, but with 


R 3 ~ 0.987 x R 1} 

R 4 = 0.176 x R 2 . 

The present calculations in Figures 1 and 2 may be compared with those 
of Hayes & Nussbaumer (1983). Our results for Ri and R 2 are up to ~15% 
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different from those of Hayes k Nussbaumer, which is principally due to the 
adoption of improved electron impact excitation rates in the present analysis. 


3. RESULTS AND DISCUSSION 


Observed values of Ri - R4 Have been measured from high resolution ul- 
traviolet spectra obtained with theTUE satellite. The objects considered (the 
planetary nebulae NGC 3918 and NGC 7662, and the symbiotic stars RR Tel 
and V1016 Cyg) are listed in Table 1, along with the IUE images used in the 
analysis. Also given axe the derived emission line ratios, which were deter- 
mined using the Goddard Regional Data Reduction Facility. We note that the 
O IV 1404.8 A feature is blended with the S IV 3i?3p 2 Vi/ 2 ~ 3s3 p 2 4 Pi / 2 line. 
However this S IV tr an sition is predicted to be <20% of the intensity of the 
S IV 3s*3p 2 P 3 / 2 - 3s3p 2 4 P s /2 component at 1406.0 A (Dufton et al. 1982). 
From the extreme weakness (or absence) of the 1406.0 A line in our spectra, 
we estimate that S IV contributes <3% to the total strength of the O IV/ S IV 
1404.8 A blend, and hence may be neglected. The spectra of RR Tel and NGC 
7662 are plotted in Figures 3 and 4, respectively, to illustrate the quality of 
the observational data, and to show the relative strengths of the O IV 1404.8 
A and S IV 1406.0 A lines. 

In Table 2 we summarise the electron densities derived from the observed 
values of Rj - R, in conjunction with the calculations in Figures 1 and 2, 
along with the adopted electron temperatures from the references listed. An 
inspection of the table reveals that for NGC 3918 and NGC 7662 the derived 
densities are compatible, with values that differ by typically <0.2 dex from 
the mean estimates ( log N e = 3.7 and 4.0 for NGC 3918 and NGC 7662, re- 
spectively). However in the case of V1016 Cyg, the electron density predicted 
from R.j is approximately a factor of 50 lower than those estimated from Ri 
and R3. This is probably due to the fact that the 1397.2 A line in the V1016 
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Cyg - spectrum is weak and shows an asymmetry in the line profile, implying 
that it may be blended. We note that the observed R» ratio would need to be 
decreased by ~20% for it to give the same electron density as Ri and R3 • 
Also listed in Table 2 are the electron densities derived from line ratios in 
species with si mi lar ionization potentials, and hence spatial distributions to 
0 IV, such as 1(4711 A)/I(4740 A) in Ar IV. The values of log N„ estimated 
from R* - R« may be seen to be very similar to those deduced from other line 
ratios, with discrepancies that average only 0.2 dex, when the R4 result for 
V1016 Cyg is excluded. In addition, we note that the observed ratios that are 
in the high density limit (Ri - R4 in RR Tel and Rj in V1016 Cyg) only differ 
from the predicted high density values by <10%. These results provide exper- 
imental support for the accuracy of the present 0 IV line ratio calculations, 
and hence the atomic data employed in their derivation. 


We would like to thank A. E. Kingston for his continued interest in the 
work. ESC is grateful to the Science and Engineering Research Council (UK) 
for financial support, while DAB acknowledges financial assistance from the 
Culham Laboratory, the Department of Physics (QUB) and the Department 
of Education for N. Ireland. AKP acknowledges partial support from NASA 
grant NAG 5-1644. The IUE archival portion of this work was supported by 
the Goddard Regional Data Analysis Facility (RDAF). This work was finan- 
cially supported by NATO travel grant 0469/87. 
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FIGURE CAPTIONS 


FIG. l.-The theoretical 0 IV emission line ratio Ri = I(2s 2 2 p 2 P 3 / 2 - 2s2p 2 
4 P x/2 )/I(2s 2 2 p 2 P 3/2 - 2s2p 2 4 P 5 / 2 ) = 1(1407.4 A)/I(1401.2 A), plotted as a 
function of logarithmic electron density (N e in cm -3 ) at electron temperatures 
of (from top to bottom), T e = 5000, 10000, 15000, 20000, 30000 and 40000 K. 


FIG. 2.-The theoretical 0 IV emission line ratio R 2 = I(2s 2 2p 2 P 3 / 2 - 2s2p 2 
4 P 3 / 2 )/I( 25 2 2 p 2 P 3/2 - 2s2p 2 4 P 5 / 2 ) = 1(1404.8 A)/I(1401.2 A), plotted as a 
function of logarithmic electron density (N e in cm -3 ) at electron temperatures 
of (from top to bottom), T e = 5000, 10000, 15000, 20000, 30000 and 40000 K. 


FIG. 3-High resolution IUE spectrum (SWP 20247) of the symbiotic star 
RR Tel in the wavelength region 1396-1408 A, where the flux is in units of 
ergcm-V 1 . The 0 IV lines at 1397.2, 1399.8, 1401.2, 1404.8 and 1407.4 A 
are clearly labelled in the figure, as are Si IV 1402.8 A and S IV 1406.0 A. 
An emission feature due to a cosmic ray event is labelled in the figure with a 
cross. 
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FIG. 4 -High resolution IUE spectrum (SWP 4106) of the planetary nebula 
NGC 7662 in the wavelength region 1396-1408 A, where the flux is in units 
of ergcm-V 1 . The 0 IV lines at 1397.2, 1399.8, 1401.2, 1404.8 and 1407.4 
A are clearly labelled in the figure, as are Si IV 1402.8 A and S IV 1406.0 A. 
An absorption feature due to a reseau mark is labelled in the figure with the 
letter R. 
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TABLE 1 


OBSERVED O IV LINE INTENSITY RATIOS 


Object 

SWP Image 
Number 

R 1 

NGC 3918 

19888 

0.23 

NGC 7662 

4106 

0.30 

RRTel 

20247 

0.15 

V1016 Cyg 

13432 

0.17 


R 2 

R 3 

R4 

0.62 

0.24 

0.16 

0.69 

0.29 

0.13 

0.55 

0.14 

0.10 

0.51 

0.19 

0.12 
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